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EMPIRICAL SPACINGS OF UNFOLDED EIGENVALUES 


KRISTINA SCHUBERT AND MARTIN VENKER 


Abstract. We study random points on the real line generated by the eigenval¬ 
ues in unitary invariant random matrix ensembles or by more general repulsive 
particle systems. As the number of points tends to infinity, we prove convergence 
of the empirical distribution of nearest neighbor spacings. We extend existing 
results for the spacing distribution in two ways. On the one hand, we believe the 
empirical distribution to be of more practical relevance than the so far considered 
expected distribution. On the other hand, we use the unfolding, a non-linear 
rescaling, which transforms the ensemble such that the density of particles is 
asymptotically constant. This allows to consider all empirical spacings, where 
previous results were restricted to a tiny fraction of the particles. Moreover, we 
prove bounds on the rates of convergence. The main ingredient for the proof, 
a strong bulk universality result for correlation functions in the unfolded setting 
including optimal rates, should be of independent interest. 


1. Introduction 

The universal behaviour of eigenvalue statistics of random matrices has attracted 
much interest over the last decades. Although random matrices have already been 
studied by Wishart }Wis28j in the 1920s, the development of the theory was pro¬ 
moted in the 1950s by Wigner |Wig58|, who used eigenvalues of random matrices to 
model spectra of complex nuclei. Montgomery’s surprising discovery |Mon73| that 
zeros of the Riemann zeta function behave statistically as eigenvalues of random 
matrices, led to a further increase of interest. In recent years, the belief has emerged 
that limit laws obtained in random matrix theory are also ubiquitous in large sys¬ 
tems of strongly correlated particles on the real line. One instance of this belief 
is the Bohigas-Giannoni-Schmit conjecture [BGS84] , stating that level spacings of 
quantum systems with classically chaotic motion should be given by random matrix 
laws. The ubiquity of certain limit laws has been established within Random Matrix 
Theory (RMT) as the universality phenomenon, which means that for large matrix 
sizes many eigenvalue statistics exhibit the same limit distributions for essentially 
different matrix models, provided these models share the same symmetry (e.g. real- 
symmetry, complex-Hermiticity etc.). These universal limits usually arise for gap 
probabilities, spacing statistics or correlations of close eigenvalues in the bulk or at 
the edge of the spectrum if the mean spacing between consecutive eigenvalues is one. 

This paper is motivated by the following problem: Assume that we consider a 
complicated (real-world) system. Based on a data set of real values obtained as a 
particular realization of that system, we want to study the appearance of universal 
RMT laws. The central question is how to detect such a universal behavior. 

From a practical point of view, the easiest and most common statistic to consider 
is the empirical nearest-neighbor spacing distribution: for an ordered iV-tuple x\ < 

■ ■ ■ < xn and an interval In we denote by <t(/jv,x) the counting measure of the 
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nearest neighbor spacings in In as 

<t(I n ,x):= 6 xj+i-xr ( L1 ) 

•Ej £-1N 


For a non-ordered x, we may define a(I^,x) := ct(In,x), where x is an ordered tuple 
build from the elements of x. The empirical nearest-neighbor spacing distribution 
of x in In is then given by the probability measure 


) x 

If / 0 °° da(Ijy,x ) = 0, then we define a(I n,x) as an arbitrary probability measure on 

R. 

Here, <j(In,x ) represents a histogram of the spacings from the data x and has as 
such a high practical relevance. Indeed, histograms of spacings obtained in numerous 
(real-world) systems have been compared to limit distributions from RMT, ranging 
from level spacings in nuclear physics (see Mehta’s classical book [Meh04] and the 
references therein) to bus waiting times in certain Mexican cities [KS031 . 

Despite its relevance, results on empirical spacings in otherwise well-understood 
random matrix ensembles are surprisingly sparse. Most results in the literature 
are available for the expected spacing distribution instead of empirical spacings. 
Briefly speaking, the expected spacing distribution is obtained by averaging over 
all realizations of ct(In,x) or d(I n,x). A more formal definition and discussion 
is given below. This preference of the expected spacing distribution is partly due 
to its direct relation to correlation functions (marginals of fixed dimension). For 
important classes of random matrices, these have particularly nice forms. This led 
to the convention of establishing local universality in terms of correlation functions. 
However, to deduce strong universality results of empirical spacing statistics, quite 
strong forms of convergence of the correlation functions are needed, e.g. uniformity 
in intervals growing with N. These requirements are often not met by standard 
formulations of universality results. 

In this paper, we prove the convergence of the empirical spacing distribution of 
unfolded eigenvalues or more general particles on the real line to a universal dis¬ 
tribution, the Gaudin distribution. The unfolding basically consists of applying 
the limiting spectral distribution function to the eigenvalues/particles. This non¬ 
linear rescaling transforms the limiting spectral measure into the uniform measure on 
[0,1] and allows for considering spacings of eigenvalues/particles from In of macro¬ 
scopic length, even the whole spectrum. Our main theorem, Theorem 12.21 states 
the uniform convergence of the distribution function of ct(In,x) towards the Gaudin 
distribution function G in mean, i.e. 




lim E n I sup 

TV-s-oo \seR 


f 


dd(lN,x) — G(s ) 


= 0 . 


( 1 . 2 ) 


We also obtain rates of convergence in terms of the length of the interval In- 

Let us define the two models of this paper. In the first model, we consider the 
eigenvalues of Hermitian matrices from unitary invariant ensembles: Given functions 
V, f : J — > R on J = [L_, L + \ n R, which is a finite or infinite interval (—oo < L_ < 
L + < oo, for precise assumptions see (GA),), we define a density on R^ by 

Pn,vj(x) := —- — JJ \xi - x j \ 2 e~ N ^= lV ^ + ^= lf ^ Xj hj(x j ). (1.3) 

Zn ’ V -J' 1 < t <J<N 


If / = 0, then we write Pn.v instead of Pn,v,o- We will slightly abuse notation 
by not using different symbols for the measure and its density. Pn,vj is the joint 
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density of the eigenvalues of a random Hermitian matrix whose distribution has a 
density proportional to exp(— Tr(NV(M) + f(M ))) with respect to the “Lebesgue 
measure” dM = nf=i dMjj rii<j<A;<jv dMjldMj k on the space of N x N Hermitian 
matrices M with the property that all eigenvalues of M lie in J. Most prominent in 
this class and arguably in the entire RMT is the Gaussian unitary ensemble (GUE), 
which is obtained by choosing V(t) = f 2 , / = 0 and J = R. 

As a second model, we will consider ensembles recently studied by Gotze and the 
second author in [ GV14] , which generalize (11.31) by allowing for different interactions 
between particles. Given Q, h : R —>• M smooth (see (GA ) 2 for our assumptions), 
we define 

P&,q{x) = -±- If Vi - ( 1 . 4 ) 

i<i<j<N 

We will call Pnq a repulsive particle system. As h is smooth, the densities Pfc q 
and Pn,vj vanish at the same order if two particles approach each other. It is 
conjectured (cf. G V1 1 ) that the universal local limit laws should only depend on 
the strength of the repulsion but be independent of the interaction at a non-zero 
distance and of the external field Q. 

Returning to spacings, let us review known results in RMT. The spacing distribu¬ 
tion has been understood best for a particular ensemble of random unitary matrices, 
the circular unitary ensemble (CUE). It was introduced by Dyson [Dys62| as a Haar 
distributed unitary matrix, which may be seen as having a uniform distribution on 
the group of N x N unitary matrices. The eigenvalues lie on the unit circle and have 
a joint density proportional to 

- e i6l \ 2 , ( 1 . 5 ) 

j<l 


where 9 \,..., On € [0, r). Note the similarity in the interaction of (11.51) and (11.31) . 
It has been shown by Katz and Sarnak IKS99| that for some limiting function G we 
have 


lim Ecue.nsup 

AT->°o sgR 


C s 27r 

/ da(I N ,(N/27r)9)-G(s) 

JO IAW| 


o (Wr 1/6+e ) 


( 1 . 6 ) 


for any e > 0, where In = [0,2itN). Soshnikov | Sos98] proved that (11.61) holds with 
a rate of 0 ((|/at|)~ 1//2+£ ) for any e > 0, where In C [0, 2irN) is such that \In\ oo 
as N —> oo, in expectation as in (11.61) and also almost surely. Furthermore, he proved 
that the fluctuations around the limit G(s) are Gaussian. 

In (fL6D . 


0(«) := E 

k> 2 


(- 1 )* 

(*:-!)! 



z j))l<i,j<k} 


Zl =odz 2 ...dz k 


is the distribution function of the so-called Gaudin distribution and S is the sine 
kernel, 


S(t — s ) 


sin(7r(t — s)) 
n(t — s) 


(1.7) 


This probability measure had already been studied before in the physics literature 
(see e.g. (Meh04j for references). In particular, the density G'(s) is given as the 
second derivative of the Fredholm determinant of the integral operator with the 
sine kernel on L 2 ((0, s)). Although this density does not seem to have a nice closed 
form expression, for many practical purposes it is sufficient to consider the so-called 
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Wigner surmise instead, i.e. P 2 (s) := ^js 2 e 4 s provides a good approximation to 
G'(s). 

Before we can further review results on the spacing distribution for invariant 
ensembles, we first need to introduce the notion of the equilibrium measure in order 
to rescale the particles. It is well-known that under very mild assumptions on V 
and J, there is a measure py on R with compact support that is the weak limit 
of the expected empirical spectral distribution of Pn,v, he. for all continuous and 
bounded g : R —> R 



Here Ejv,y denotes the expectation w.r.t. Pn,v- The measure py is the unique 
probability measure among all Borel measures on J which minimizes the functional 


p i-A 



1 dp(t)dp(s) + 


V(t)dp(t). 


Throughout our paper we will assume that V is convex and real-analytic on a neigh¬ 
borhood of the support of py (see (GA)i). Then py has a density, which we will 
abusing notation also denote by py, strictly positive in the interior of the support 
of py. 

In invariant ensembles the results for the spacing distribution are by far weaker 
than for circular ensembles. For instance, until the recent [Schl5| by the first author, 
only the absolutely continuous intensity measure E<t(/jv,x) had been considered for 
these ensembles. We recall that if x is random, then <j(In,x) is a random measure 
and its intensity measure is defined as 


IB 


dE<r(/jV ) x ) := E / da(I]y,x ) 


(1.9) 


IB 


for any measurable set B and E denotes expectation w.r.t. the probability measure 
underlying the random variable x. We will call this measure the expected spacing 
measure. 

To our knowledge, the first rigorous result on the spacing distribution for unitary 
invariant matrix ensembles is due to Deift et al. DKM + 99] . Let a E supp (py)° and 
tjy > 0 such that tjy —> oo and tjy/N -A 0 as N -A oo. Setting In = [a — tjy, a + f/v], 
it has been shown in |DKM + 99] following the method of Katz and Sarnak that for 
real-analytic V and / = 0 in (|1.3I) . we have 


lirn 1 ( dE N ya(l N , Np v (a)x) = G(s). (1.10) 

N ^oo \I N \pv(a) J 0 

We observe that (11.101) clearly shows universality, as the r.h.s. of (11.101) does not 
depend on V. Let us also remark that (11.101) expresses the weak convergence of the 
distribution function of the asymptotically normalized expected spacing distribution 
to the Gaudin distribution. Furthermore, due to the continuity of G, (II . 10p actually 
holds uniformly in s. 

An analogous result for certain Hermitian Wigner matrices was proved by Jo¬ 
hansson in [•Ioh01| . Universality was proved for large classes of Wigner matrices 
and invariant ensembles in terms of the expected spacing measures by Bourgade, 
Erdos, Schlein, Yau, Yin et al. (see BFV1 1 and the references therein). In par¬ 
ticular, they show vague convergence of the asymptotically normalized expected 
spacing measure. Moreover, the limiting distribution has been proved to be univer¬ 
sal for large classes of real-symmetric, Hermitian or quaternionic self-dual random 
matrices as well as for general /3-ensembles, which are variants of (11.31) in which 
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a parameter j3 > 0 replaces the exponent 2. The universal limit depends on the 
symmetry type or on /3, respectively. 

Recently, there has been quite some interest in the distribution of a single spacing. 
This was initiated by Tao, who proved in [Taol3 j that the Gaudin distribution is 
also the limit in law of a single spacing in the bulk of the spectrum, say of ®[tv/2J+i — 
x I tv/ 2 j • This shows that for expected spacings, an average over an increasing number 
of spacings as in (11.101) is not necessary to obtain the limiting distribution. Tao 
proved this for Hermitian Wigner matrices and the result was later extended in 
[EY12j and [BFG15 ] to all symmetry classes and /3-ensembles. Moreover, the results 
in jEY12j and |BFG15j allow to consider statistics of the form (for ordered Xi) 

Ejv Y] 7| Xi)), (1.11) 

I I iel 


where /M/.j '■= Hv(Qi) and qi is the z/lV-quantile of the distribution jiy. The test 
function g is assumed to be smooth and of compact support, thus determining 
vague convergence. The index set / is assumed to contain only bulk eigenvalues in 
case of |EY12| and is arbitrary in [BFG15] , Both works show that in the large N 
limit, (11.111) become independent of V, thereby showing universality. However, to 
deduce convergence to the Gaudin distribution via [ Taol3] , the set / has to fulfill 
eN < min / < max / < (1 — e)N for some e > 0. 

It is instructive to compare the statements (11.111) and (11.21) . In (11.111) as well as 
in (11.10ft . by first taking the expectation, the spacing statistics are averaged over 
all realizations, thus making it non-empirical. After that, this average is compared 
to the limiting quantity. In (11.21) . the empirical statistic is first compared to the 
limit. The error of that comparison is then shown to vanish in L 1 . It is further 
important to note that for statements like (11.21) . the number of eigenvalues in the 
statistic necessarily has to increase with N. Indeed, the smooth function G can only 
be approximated well by the step function f Q da(I^r,x) if the number of steps goes 
to infinity. The convergence of (11.111) for a finite set I is only possible since taking 
the expectation means averaging over infinitely many spacings (from all different 
realizations). 

Let us turn to empirical spacings. In |Schl5| the first author of this work shows 
under certain assumptions on the Christoffel-Darboux kernel (see (12.181) for a defi¬ 
nition of the kernel) that for unitary invariant ensembles 


lim Eyv v sup 

AT-S-OO ’ seR 


1 


/o w(.a)\I N \ 


da(Ijy,Ng,y(a)x) — G(s ) 


= 0 


( 1 . 12 ) 


with /tv, a as in (11.101) . This result is quite analogous to (11.61) for the CUE and was 
also shown in [Schl5j for eigenvalues of real-symmetric and quaternionic self-dual 
ensembles. 

However, this result has the drawback that the considered spacing statistics are 
very inefficient in an empirical sense in that only a tiny fraction of the eigenval¬ 
ues is used for the statistics. Indeed, the expected number of rescaled eigenvalues 
N/iy(a)xj in /tv = [a — t^, a+tjy] is about 2Uv and thus by the condition tjy/N —>• 0, 
the fraction of used eigenvalues even goes to 0, as N gets large. This unsatisfying 
situation is due to the scaling of the eigenvalues 


N fiy(a)xj, j = 1,..., N, (1.13) 

which we will call localized around the point a. To obtain a universal limit, the 
eigenvalues have to be rescaled such that their mean spacing is (asymptotically) 
one. This is achieved by the rescaling (11.131) in a small vicinity of a point a in the 
bulk of the spectrum, where the mean spacing is (/ry(a)A^)" 1 + o(N~ 1 ). But as the 
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equilibrium density /ay (a) is generically not constant in a (in contrary to the CUE), 
a linear rescaling such as (11.1311 does not allow to consider spacings of eigenvalues 
which are spread over an interval of macroscopic size. To overcome this problem, 
we will in this article unfold the eigenvalues with the distribution function of fly. 
This natural non-linear rescaling produces an ensemble with constant equilibrium 
density. 


2. Statement of Results 

Before we give a precise formulation of our results, we specify our assumptions for 
the ensembles (11.31) and (| 1 .41) . Except for the support J, the class of ensembles Pfc q 
formally subsumes the unitary invariant ensembles Pn,v- However, our assumptions 
depend on the specific type of ensemble and therefore we will consider two different 
sets of assumptions. For Pn,v,j these assumptions are 

(1) V, f : J M are real analytic, J = [L_,L+] n M with —oo < 
L_ < L + < oo. 

(2) inf te j V"(t) > 0. (strict convexity) 

(3) limui^oj, V(x) — cf(x ) = oo for some c > 0 in case that J is 

(GA) j unbounded. 

(4) L_, L + do not belong to the support of the equilibrium measure 
Mv- 

If / = 0, then (2) can be replaced by the weaker condition 
(2’) V' is strictly increasing. 

For the ensembles Pfo q, we make the following assumptions. 


(1) Q, h : M —>• R are real-analytic and symmetric around 0. 

(GA) 2 ( 2 ) a Q := infteK Q"(t) > 0. 

(3) h is a Schwartz function with exponentially fast decaying Fourier 
transform. 

Under (GA) 2 , it has been shown in [GV14j that for given h, for each Q with 
aQ large enough (depending on h), there exists a probability measure fiQ such that 
(11.81) holds with //q replacing fly and IEreplacing E jv,v- Furthermore, /Xq has 
also the generic form (14.471) for V := Q + h * /Xq, where * denotes convolution. 


Remark 2.1 (On the assumptions). One could without substantial changes also intro¬ 
duce the external held / to the model Pfc q ■ For the sake of notational convenience 
this is not done here. The model could also be studied on an interval J, which might 
not be the whole line. This would require a condition analogous to (GA)i (4) for 
Pj^r q. For repulsive particle systems the equations (I4.42|) , which determine the end¬ 
points of the equilibrium measure, depend on the equilibrium measure itself, which 
makes checking an analog of (GA)i (4) more complicated than for invariant ensem¬ 
bles. Nevertheless, we will show in our analysis that a truncation to large enough J 
only produces an asymptotically negligible error. 
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To define the scaling of the particles, for Pn,vj let Fy denote the distribution 
function of the equilibrium measure py and consider the unfolded eigenvalue 

Xi := NF v (xi). 

If x is a random configuration sampled from Pn.vj- then x is a point process on 
[0, N] with asymptotically constant density and mean spacing 1, at least in the bulk. 
If x is distributed according to the repulsive particle system Pjf q , let Fq denote the 
distribution function of jj h q and define 


Xi ■= NF%( Xi ). 

Now, our main theorem on the spacing distribution reads as follows. 


Theorem 2.2. Let Ijy C [0, IV] be a sequence of intervals with 

liminf — : dist(/jv, {0, N}) > 0 and lim |Ijv| —» oo. 

TV—>• OO TV N — ^OO 


(1) Let Pr\r be either Pn,vj satisfying (GAh or Pjf q satisfying (GA ) 2 with 
cxq large enough (this depends on h only) and h negative-definite. Then for 
any e > 0 


Etv ( sup 
\seR 


da (In, x) — G(s) 



(2.14) 


(2) If Pn = Pnq w 'dh h not necessarily negative-definite, then for a.Q large 
enough ([2.141) holds with the O-term being replaced by o(l). 

In (1) and (2), (t(In,x) can be replaced by |Uv| -1,J (-^Ar, x) without altering the result. 


Remark 2.3. 

(1) In the language of mathematical statistics, (I2.14|) implies that a(I n,x) is an 
asymptotically consistent estimator for the Gaudin distribution, considered 
in the Kolmogorov metric. Note that similar results on the expected spacing 
distribution like (11.101) only show the asymptotic unbiasedness. For applica¬ 
tions, it would be favorable to have the unfolding with F := Fy or F : = Fq 
being replaced by an unfolding based solely on the empirical values. We 
remark that using the empirical distribution function as a naive estimator 
for F would result in a non-random x. However, we expect that a smoothed 
empirical distribution function should yield the desired. In fact, this is a 
typical procedure in statistical analysis. It will be considered in a future 
work. 

(2) To our knowledge, Theorem 12.21 is the first result on rates of convergence 
to the Gaudin distribution for invariant ensembles. For the simpler CUE a 
rate of 0(\In\~ 1 ^ +£ ) for any e > 0 was derived in [KS99] resp. a rate of 
0(|/jv| -1 / 2+e ) was shown in [S os98] . Numerical experiments (cf. [KS13] ) 
suggest that the optimal rate for the GUE is |/at| - 1 / 2 , possibly with some 
logarithmic factor. The dependence of the rate on \In\ reflects the fact that 
necessarily a growing number of empirical spacings has to be considered in 
order to obtain convergence. 

(3) For negative-definite h, an exact representation of P), q in terms of determi- 
nantal ensembles will be derived in Section [5j which allows to transfer rates 
of convergence from the unitary invariant ensembles to the repulsive particle 
systems. For general h, only convergence can be shown, see Remark 15.21 for 
more details. On the other hand, if h is positive-definite, then it suffices to 
have > sup tgR —h"(t). 
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Abandoning any rate of convergence, we an also deduce the result of Theorem 12.21 
for I N = [0, N\. 


Corollary 2.4. With Pjv as in Theorem 1 2. 21 (1) or (2), we have 


lim E^r sup 
N ~> °o VseR 


cZ<t([0, N], x) — G(s) 


= 0. 


(2.15) 


Remark 2.5. Note that in (|2.15l) . edge spacings are included. Although correlations 
between eigenvalues at the edge are not given by the sine kernel, the number of edge 
spacings is relatively small and thus does not change the limit distribution. 


The next corollary shows a much better rate of convergence for the expected 
spacing distribution. We believe that this rate is almost optimal. 


Corollary 2.6. Let Pn,In be as in Theorem, 1 2. 21 (1). Then for any e > 0 


sup 

sSR 



1 

\ j n\ 


dE]\rcr(I^r,x) — G(s) 


o (i/vr 1+£ ). 


Remark 2.7. Recall from the discussion around (11.111) that a similar result should 
also hold for intervals with |ijv| of order 1, probably with rate 0{N~ l+£ ). As 
Corollary 12.61 is merely a byproduct of an analysis which necessarily deals with 
growing intervals, we will not pursue this here. 


A major ingredient to the proof of Theorem 12.21 is a strong form of bulk univer¬ 
sality for correlation functions, which should be of independent interest. To state it, 
let us recall the notion of correlation functions. 

For a probability measure Pj\r(x)dx on R^, invariant under permutations of its 
arguments, the fc-th correlation function is the map /?/(, : R fc — >• R, 

i, n\ r 

.. ,t k ) := ——— Pisr(ti,...,tN)dt k+1 ...dtN. (2.16) 

{N-k)\J R N-k 

Note that in previous works of the second author, the correlation functions are 
defined as the /s-dimensional marginal densities of Pn and therefore differ from 
(|2.16l) by the factor N\/(N — k)\. The definition in (|2.16l) is more convenient for 
dealing with sums and will be used throughout this article. 

A crucial fact for many computations and universality proofs is that unitary 
invariant ensembles are determinantal, that is 

(2-17) 

such that the analysis boils down to studying the so-called Christoffel-Darboux ker¬ 
nel 

N—l 

K N yj(t,s ) := 2Pi.^WPi,Jv(s)e" JV/2(V(t)+VW)+1/2(/(t)+/W , (2.18) 

3=0 

where Pj t N,j = 0,1,... are the orthonormal polynomials with positive leading coef¬ 
ficients to the weight e ~^^^ NV on J. 

To formulate a certain uniformity in the field /, let for a domain D C C, (Xd, || • 
||£>) denote the Banach space of functions / : D —>• C which are analytic, bounded 
and real-valued on D DR. Here || • ||£> denotes the sup-norm on D. 


Theorem 2.8. Let In C [0, IV] be such that lim inf n- 
(1) Let V, f satisfy (GA)\. Then 


, dist(/;v, {0, N})/N > 0. 


1 


-K 


N,V,f 


V 




-1 

V 


Np V (Fy\jj)) 

where the error term is uniform for t,s € In- 


sin(7r(f — s)) 
ir(t — s ) 


+ o 
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( 2 ) 


(3) 


Assume in addition to (1) that J is compact and let 0 < 77 < 1. Then there 
is a complex domain D D J such that (1) holds uniformly for f £ Xp with 
||/||d < N v if 0(1/N) in (1) is replaced by 0(N 71 ~ 1 ). 

Let B/fj be the k-th correlation function of Pn with Pn as in Theorem \2.2\ 
(1). Then, abbreviating tj := Fy l (tj/N) in the unitary invariant case and 
tj := (FQ)~ l (tj/N) in the case Pn = P^q, as well as writing p for py and 

T h Q 


respectively, we have for any e > 0 


1 


N k UU . 


—R%(ti, ... ,4) = det 


sin(7r(tj - tj)) 
7 t(U - tj) 




+ O (N~ l+£ ) 


with the O term being uniform for t\,... ,t k € In- If Pn = Pn,vj , then the 
statement is valid for e = 0. For Pn as in Theorem \2.2\ (2), the statement 
is valid with the O term replaced by o(l). 

Remark 2.9. The convergence of the Christoffel-Darboux kernel of some determi- 
nantal ensemble to the sine kernel is a very typical result in RMT and the content 
of numerous papers in the field. We will mention here only very few seminal papers 
and refer to [Kuill] for an overview instead. A first universality proof was given by 
Pastur and S hcherbina jPS97] (see also [ PS08} ) for sufficiently smooth V. Deift et 
al. DK M + 9 9j showed universality for real-analytic V using Riemann-Hilbert tech¬ 


niques and more recently Levin and Lubinsky |LL08| established it under very mild 
assumptions on V using complex analysis. 

In the existing literature, the kernel is usually considered in the localized scaling, 
that is Kn(o, + N *( a \ , a + N ^( a ) ), where a is in the bulk of the spectrum and p is 
the limiting spectral density. Then t and s are typically assumed to lie in some 
fixed compact set, that is their distance is bounded in N. In the recent [KSSV14] , 
convergence is shown under the assumption \t — s| = o(N) with the rate 0(( 1 + \t\ + 
|s|)/iV), which is optimal in the localized scaling. To our knowledge, Theorem 12.81 
is the first version of bulk universality for unitary invariant ensembles which does 
not require all eigenvalues to lie in a vicinity of some point a. Moreover, the rate in 
part (1) of the theorem is optimal, as can be seen for instance in (14.491) . 

Theorem 12.21 will be deduced from the following more general result. 

Theorem 2.10. Let for each N, In C [0, 00 ) be an interval and Pn(x)Ix be a 
probability measure on M. N , invariant under permutation of the coordinates. Let 
R k N denote the k-th correlation function of PN(x)dx, defined in (12.161) . Further let 
Cq > 0 denote a constant such that the following conditions are satisfied: 

(1) For all N € N we have 


sup 


J L-R%(t 1 /N,...,t k /N) 


<c k . 


(2) There exists kn > 0 with lmiAr_ 


sup 

£l,...,tfc€/jV 


j^R k N (t i/ N , 


,t k /N)-S k (t 


kn = 00 such that for all N € N we have 

1 

k n 


i) 


i Ik) 


< k k/2+1 • C 0 fc -, 


where for S as in (E21). S k is given by 

R k (l i) ■ ■ ■ ■ det [<S(tj tj)]|<i j< k ■ 

Then for any e > 0 

Ejv (sup f da(lN, Nx) — G(s) ) = O (mm(\lN\, kn) 

VseR Jo J ^ 

The same result holds with ct(In,Nx) being replaced by \In\~ 1 o'(In,Nx)- 































EMPIRICAL SPACINGS OF UNFOLDED EIGENVALUES 


10 


The preceding theorem can also be used to show universality of the spacing dis¬ 
tribution in the localized scaling. 


Corollary 2.11. Let Pm be as in Theorem \2.2\ (1). Set pL = fj,y if Pn = Pn,vj and 
H = Hq if Pn = Pm q ■ Let a be such that n(a) > 0 and Im C J be a symmetric 
interval with center a and linpv-^ool-T/vl = oo but linuv_ > . 00 |/Ar|/A’ = 0. Then 


Ejy ( sup 


/ 

Jo 


da (In, N/j,(a)x) — G(s) 


= O (mind/jvl, N/\I n \Y 


If Pn is as in Theorem \2.2\ (2), the last statement is valid with the O term being 
replaced by o(l). In either case, the statements remain valid if d(I n, NpL(a)x) is 
replaced by (|/jv|M(a)) _1 <7(/jv, Npt(a)x). 


Remark 2.12. Note that for \Im\ VN, Theorem 12.21 provides a better rate of 
convergence for the unfolded particles than Corollarv l2.11l for the localized particles. 


We will finish this section with some concluding remarks. The repulsive particle 
systems Pjf q appeared first in a more general setting with many-body interactions 
in |BdMPS95] . where under some convexity condition on the additional interaction 
results on the equilibrium measure were announced. Further results associated with 
global asymptotics in such models can be found in [BGK13] and [CGZ14j . Local 
bulk universality was proved for P^ q hi [ GV14 ] and for the /3 variants in [Venl3j . 
Edge universality and fine asymptotics of the largest particle have been considered 
for Pfc q in the recent [KV15| . 


The paper is organized as follows: After a brief outline, the proof of Theorem 
12.101 will be given in Section [31 Here we follow the method developed by Katz and 
Sarnak in [I KS 99] , streamlining and optimizing in order to obtain better rates of 
convergence. Theorem 12.81 (1) and (2) will be proved in Sectional The proof of 
Theorem 12.81 (3) for the repulsive particle systems relies on a non-trivial reduction 
to the unitary invariant case. An outline of the method from j [GVl4] and the proof 
of Theorem 12.81 (3) are contained in Section [5j The proofs of Theorem 12.21 and the 
corollaries are given in Section [6j 


3. Investigation of the Spacing Distribution - Proof of Theorem 12.101 


Theorem 12.101 will be proved first with the asymptotic and non-random normal¬ 
ization |ijv|, the case of the exact but random normalization will be discussed at the 
end of the proof. Moreover, let us note that in this case the statement of Theorem 
12.101 is trivial if |/jv| is bounded in N. Hence, we will from now on assume that 
\In\ —>• oo, as A-> oo. 

Let us furthermore make some notational remarks. By C we denote absolute 
positive constants that may change from line to line. Finally, note that we often 
suppress certain IV-dependencies. 

A major disadvantage of g(Im,x ) is its dependence on the ordering x\ < ■ ■ ■ < 
xn, which prevents an efficient use of correlation functions. This problem can be 
circumvented by using the measures 


7 k (lN,x ) 


— y 

|/Ar| u< ~<i k , 


^(maxi< j 


<k 


—mini<j<i, Xij), 


which are symmetric and fulfil the relations 
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T^-T r da(I N ,x) = V(-l) fc [ S d 7 k (I N ,x), N £ N, 

I-‘AT | Jo k=2 JO 

(-l) m J-r [ S da(I N ,x)<(-irf2(-l) k [ S d^ k (I N ,x), m < N. 
l-'ivl Jo k=2 Jo 


(3.19) 

(3.20) 


The proof of Theorem 12.101 consists of three steps. The first step establishes the 
point-wise convergence 


lim Ejv ( [ Tv— rda(I]y, Nx)) = G(s) 
A->-oo Wo / 


(3.21) 


and bounds the difference of j^\da(I]y, Nx) and G(s) in terms of the variances 

of the 7 fc ’s. In the second step, these variances are estimated. From this we can 
deduce a bound on 


Eat 


I-^atI 


f da(lN,Nx) — G(s) 
Jo 


(3.22) 


(see Corollary 13.411 . 

The difference between 13.221) and the quantity to be estimated in Theorem 12.101 
is that we need to take the supremum over all s before taking the expectation. This 
issue is addressed by considering (13.221) at a (growing) number of nodes s, rather 
than at a single s. The respective results are provided in Section [3731 


Before we turn to the proof of Theorem 12.101 we note an important estimate for 
power sums. We will frequently encounter sums of the form )T)fc =2 a kZ k with L and 
z growing in N and with different sequences (a*,). To provide a unified and efficient 
treatment of these sums, the following simple lemma will prove useful. 

For an entire function /, recall that / is said to be of finite order if the inequality 

max|/(z)| < e rK 

\z\<r 


holds for all r large enough and some k < oo. The greatest lower bound of all such 
k is called order of /. If / has a power series expansion Y1T=o a k zk > then the order 
p can be found via 


p = lim sup 

k —^oo 


k log k 
log|o fc |' 


(3.23) 


If / has finite order p, then it is said to be of finite type, if 

max|/(z)| < e^ rP 
\z\<r 


holds for r large enough and some finite The greatest lower bound v of all such 
C is called the type of / and can be determined via 

( vep) x l p = limsupfc^lafcl 1 ^. (3.24) 

k—¥ oo 


Lemma 3.1. Let f(z ) = cik\z k be an entire function of order at most 2 and 

finite type. Let (Ln)n, (8n)n, (^Jn)n be sequences with Ljv > 0, Ljsr = o(logMjv) 
and 0 < 5n < y/fogMfj. Then 


k =0 


om) 


for any e > 0. 
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Proof. Let e > 0 and v denote the type of /. For any K > 1, we have the trivial 
estimate 


k =o 


k =o 


< K Ln J>fcl (%) < K LN f(5 N /K) < CK LN e u & < CK L *Mf. 


K 


Choosing K = K(s) large enough, the lemma is proved. Here we used that K Ln is 
of subpolynomial growth in Mjy, due to our assumption on Ln- □ 


3.1. The convergence of Eat ^ fy Nx)J . We turn to the proof of (13.211) . 

For s > 0 and t\,...,t k , we denote by Xs,i N the function 

Xs,l N (h,...,t k ) ■= l(o, s ) ( max ti — min U ) TT ti N (ti). 

\i=i,...,k i=l,...,k ) ^ 

To select certain entries of a vector x = (x ±,..., xn) € M. N we use the notation 
xt ■= (x h ,.. .,x ik ), T = {ii,... ,4}, 1 < h < ... < i k < N. 

With this notation we can rewrite 


J q d~/ k (I N ,Nx) = j^-j- Xs,i n (Nx t ) 


(3.25) 


Tc{l,...,Af},|T|=fc 


and we obtain 


( I f d"y k (I N ,Nx) 


J Xs,i N (h, ■ ■ ■, t k )R$f(ti/N,..., t k /N)dt\ ... 


\I N \k\N k 


dt k . 


(3.26) 


The following lemma establishes the convergence of the terms Eat (fy d 7 fc (-Tv, -(Vx)) 
and further provides a useful estimate for the proof of Theorem 12.101 The proof is 
essentially contained in [Schl5j and we revisit the arguments in order to adjust them 
to the current setting. 


Lemma 3.2. Let the conditions of Theorem, \2.1 fh he satisfied. 

(1) For k >2 we have 

Ew ( [ d^ k (I N ,Nx)) = [ S k (0,z 2 ,...,z k )dz 2 ...dz k 

\J 0 / Jo<z 2 <...<z k <s 


1 


(fc-1)! \\I N \ 


1 


+ tt - tt rS k O ( —— ] Cy + s fc_i Ci 


-i r,k k k/2+1 1 

0 [k- 1 )! k n ’ 


(3.27) 


where the O term is uniform for s € R and k € N. 

(2) For An := min( \In\- kn), 0 < s := sag sn = 0{\J log An) and L := Ln € N 
such that ^/log An = o{Ln ) and Ln = o(logHAr), we have for any e > 0 

G(s) - j^—\ J da(I N ,Nx ) <£E,(/ d^ (I n , N x)^ - J d^ k (I N ,Nx) 
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Proof. In order to prove (1), we consider (13.261) and use the uniform convergence of 
the correlation functions in (2) of Theorem l2.101 i.e. we use that -^R k N (h/N,...,t k /N) 
S k (ti,..., tk) + k k / 2+l ■ Cq ■ -^7 uniformly on In- We further use the obvious estimate 


~ W=P- (3 ' ; 

to obtain 

E„ ([ N*)) - ^ l x,,i„(t)Sk(t)dt + 


The translational invariance of S k and the change of variables z\ = t\,Zi = tj — t\, i = 
2 ,... ,k lead to 


\I N \k) 


f 

J R* 


Xs,i N (t)S k (t)dt = 


l-Twl 


f 0^- z 2 —••• — z fc — s 

2 1 6/ JV- z l+ z j^ I N>3= 2 f 


S k (0,z 2 , ■ ■ .,z k )dzi ... dz k . 


Observe that in the latter integral, we integrate over z\ from In except for an 
interval which has at most length s. The estimate in (1) is then obvious from 
sup teI k |Sfc(t)| < Cq (see Theorem 12.101 ('ll). Observe that statement (1) together 
with (13.201) already implies (13.211) . 

To show (2), we first introduce the notation 


E(s, k) := / Sk(0,z 2 ,...,Zk)dz2...dz k . 

J 0<Z2<---<Zk<S 


(3.29) 


The idea is to use (13.201) and bound G(s) and J Q S da(lN, Nx) from above and from 
below by alternating sums over E(s,k) and dr/ k (lN, A T x), respectively. Then we 
obtain for L € N 


G(s)~ 


1 


l-^vl 


da (In, Nx) 


< Y QT d^ k (I N , Nx)^j - J* d 7 k (I N , Nx) 


+ Y En dry k (I n , Nx)^) - E(s,k) 


+ E(s, L) + E(s, L + 1). 


Introducing the notation s := max(l,s), we conclude from ()3.281) 

fiL +1 pL 

E(s,L) + E(s,L + l)<jf—^- l 
and using (1), we arrive at 

G(s) - — J da(I N , Nx) <Y E ^(^J d^^N, Nx)^j - J d^ k (I N ,Nx) 


+ 


Cn +l S L 


+ o 




E 


k=2 
S k Cn 


+ — Y,* k ~ 1 Co 

K AT U 


(L~ 1)! \\I N \) ^ 2 (k-l)\ K Nk=2 


fck/ 2+1 

W 1 w 


CjL+1-L 

Now, under our assumptions on the growth of s and L, °(l/ J ^A r )- The sum 


f „krdk 


L—l 


k =2 


(k-l)\ 


sc 0 y : ° 


k =1 


k\ 
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is 0{A £ n ) for any e > 0 by Lemma 13.11 applied with f(z) = e c ° z . To deal with the 
remaining sum, first observe that the series 


OO 


/(*):= 

k =0 


j^k/2+2 

kl 


converges absolutely and defines an entire function. Using (13.2311 and Stirling’s 
formula, we readily compute its order as p = 2 and using (13.2411 and again Stirling’s 
formula, its type as v = C^e/ 2. Thus Lemma 13. 1 1 finishes the proof. □ 


3.2. The variance of / d^y k (Iisr,Nx). Taking expectations in Lemma (3~21 (2 ). we 

Jo 

arrive at a sum of expected absolute differences of drf k {Ij\i, Nx ) to its expectation. 
We will estimate this in terms of the squareroot of the variance of fg d^ k {lNi Nx), 
which we now bound. 


Lemma 3.3. In the situation of Lemma \3.2l there exists a positive constant C such 
that for k € N, k > 2 and N sufficiently large we have 


Var 



^ C(2s) 2k Cl k 2 k 
~ I^l(fc-l)! 


+ 


( 2 k) k+1 

{k\y 


Cl k k 2 0 



s 2k 


Proof. In order to calculate the variance of fg d"y k (lN, Nx), we first consider the sec¬ 
ond moment Ejv ((Jq dr/ k (lN, Nx))~^J using the representation (13.2511 . We expand 


E TV 


rs \ 2' 

d^ k (I N ,Nx) 


’0 


2k 

^E E L Xs,i N (Nt T )Xs,i N (Nt M )RN( t ) dt - ( 3 - 30 ) 

I N ' l=k 4M jv 


l=k T,MC{1,...,JV} 

\T\=\M\=k,TUM=l 


First, we consider the inner sum in (13.3011 for l = 2k, i.e. in the case that T, M C 
{1,...,N} satisfy \T\ = \M\ = k and TnM = 0. Since there are = 

fc! 2 (jv- 2 fc)! suc ^ se ^ s ) we obtain by the symmetry of the correlation functions 

TT~]2 \ Xs,i N (Nt T )xs,i N {Nt M )RN{t)dt 


\T\ = \M\=k,TnM=(D 


1 

Jkiy 


lT 1 2 AT2k I Xs,I N (h, ■ ■ ■ ,tk)Xs,I N (t k+ 1, • • • ,t 2 k)R 2 N (tl/N, .. .,t 2k /N)dti ... dt k . 
\1n\JX jR2k 

(3.31) 


We consider the terms with l < 2k in (I3.30H later and observe that by (13.2611 with 
t! := (ii,.. .,t k ),t" := (4+i,.. .,t 2k ) we have 

(e (jT d~f k (I N ,Nx)^ 

= (fc ; )2 J R2k Xs,iAt')Xs,iAnR k N(t'/N)R k N (t"/N)dt 1 ... dtk. ( 3 -32) 

To calculate the difference of (13.3111 and (|3.32H we use 


-^{R 2k {t'/N,t"/N) - R k N (t'/N)R k N (t"/N)) 

= S 2k (t',t") - S k {t')Sk{t”) + (2 k) k+1 Cl k O (—) . 

\KN J 
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By (13.281) we obtain 

( fc l )2 \ lN \l N 2k f R2k Xs,i N {t') Xs ,i N {f){R^{i! /N,t" /N) - R k N (t' /N)R k N {t" /N))dt'dt n 

= 77712 TT12 [ ~ S k {t')S k (t"))dt'dt" 

(k\y \i N \ z .y R 2 fe 


+ 


( 2*0 


fe+i 


-Cn k k 2 0 ( — ) s 2k 


1 

k n 


(k\) 2 0 

We now claim that 

S 2 k(t', t") - S k (t')S k (t") < 0, t', t" g M fc . (3.33) 

If two components of (t' ,t") are equal, then S 2 h{t',t ") = 0 and the claim is trivially 
true, as S k > 0. If all components are distinct, then 

'sin(7r(t n - t m )Y 


t m ) ) i <n,m<j 

is a positive-definite matrix (its principal minors are exactly Si, S 2 , ■ ■ ■, Sj-i > 0). 
Now, (13.331) follows from Fischer’s inequality. With (13.331) we can further estimate 

(fc j )2 m l N 2k J R2k Xs,i N (t')Xs,i N (t")(RN(t’/N, t"/N) - R k N (t'/N)R k N {t"/N))dt'dt” 


< 


( 2*0 


fc+i 


C% k k 2 0 ( — ) s 


(k\) 2 0 

So far, we showed 


1 


k n 


—2k 


Var f f d-f k {I N ,Nx)] < — ^ V V [ X.sJ N {Nt T )Xs,i N (Nt M )R^(t)dt 

W° / \ lN \ T,Mc7t..,iV } J ^ N 


\T\ = \M\=k,TUM=l 

(2 k) k+1 3k 2 /J_\ 2k 

+ (k\) 2 C ° k0 {K N ) S • 

We continue to consider the integrals in the double sum. For l € {k ,..., 2k — 1} 
and sets T, M C {1,... ,7V} with |T| = \M\ = k and T U M = I (i.e. T and M have 
a non-empty intersection) we have (using the symmetry of R k N and jttR 1 n < C l 0 ) 

[ Xs,I N (Nt T )Xs,I N {NtM)RN(t) d t 
J R N 

(N — l)\ 1 

IV! 


< 


jji / l(o, 2 s) Qrnax^tj “ .min^ R l N (ti/N,.. .,t t /N)dti ...dt t 


< 


l\I N \{2s) 2k C 2k 


N\ 


(3.34) 


We observe that for given l there are Cjv,z,fc := (^) ([.) ( 2k -i) se ^ s T,M C {1,... ,N} 
with |T| = \M\ = k and T U M = l by some easy combinatorial argument. Hence, 
by ( | 3.34 | ) and -^(N - l)\C N ^ k = ( 2k -i)l(i-k)\* ’ we obtain 

j 2 fc —1 . 

rrni E E L Xs,I N (tT)Xs,I N (tM)RN,N{t/N)dt 

' l= k 

\T\ = \M\=k,TUM=l 
1 2/c-l 

i T (2S) 2t C 0 2 ‘ £ 


< 


|lv| 




(2 k - l)\{l - k)\ 2 ' 
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By the easy calculation 
2k-i z 


§ (2 k-l)\(l-k)\ 2 k\^ l\ \l 


^ k —1 

ctE 


l + k fk 


< 


2k 

T\ 


k -1 

E 

1=0 


< 


we get 


^ 2fc—l .. 

I^nI 2 ^ -UK'’ 


Xs,i N {t T )Xsj N (t M )RN(t/N)dt < 


Ij - (k- 1)!’ 
{2s) 2k Cl k 2 k 


k T,MC{1,...,N} 

\T\ = \M\=k,TUM=l 


M (k — 1)! ’ 


Summarizing, we have shown 


Var dj k (lN,Nx) 


< -^—(2s) 2k C$ k 2 k 1 + Cl k k 2 0 


l J v| 


(k - 1)! (k\) 2 


1 

At W 


s 2fc . 


□ 


/From Lemma 13.31 we can already derive an estimate on the expected deviation 
of the spacing distribution from its limit at a given point s. 

Corollary 3.4. Let the conditions of Lemma \3.°2\ be satisfied. Then we have for any 
£ > 0, s > 0 




1 


\ j n\ 


f da(lN, Nx) — G(s] 
Jo 


= O 


1 


1 / 2 —£ 


A 


N 


Proof. By statement (2) in Lemma 13.21 it suffices to bound 

Ejv d^ k (I N ,Nx)^j - J° d 7 k (I N , Nx)) IpE J Var (f d'yb^N, Nx) 

\k=2 / k=2 V 

The subadditivity of the square root together with Lemma 13.31 give 

t = o (-4=) 

Similar to the proof of Lemma 1 3.2 [ we can apply Lemma 13. II with the functions 

OO 

aw — E 


k=0 

oo 


Mz) = ^ (2 k)^ l)/2 C k c k/2 kz ^ 


k =0 


k\ 


which are both of order 2 and finite type, as can be checked easily using (13.231) and 

4322 ). ' * □ 

3.3. Completing the proof of Theorem 12.101 The idea for the rest of the proof 
of Theorem 12.101 is to replace the supremum over s with a maximum over a finite 
set of certain nodes s t . Then, we can choose the number of these nodes, growing 
with N in such a way, that the error estimates lead to the error claimed in Theorem 

EM 


































EMPIRICAL SPACINGS OF UNFOLDED EIGENVALUES 


17 


Let M = M(N ) € N and let Sj, i = 0,..., M, denote the i/M-quantile of G, that 


is 


^(si) = —, i = 0,...,M. 

The existence of these nodes is ensured as G is continuous by definition, increasing 
and lim^oo G(t) = 1. Now, set 

A(s,I n ,Nx) '■= — ( da(I N , Nx) — G(s), A m (In,Nx):= max \A(si,I N ,Nx)\. 

Furthermore, denote the largest non-trivial node by 

6 n := max(l, s M -i)- 

It is known (cf. [KS99I Proposition 3.1.9]) that the Gaudin distribution has sub- 
Gaussian tails, that is, 

1 - G(s) < Ae~ Bs2 (3.35) 

for some A, B > 0. This implies that Sn fulfils 

S N = O(VlogM). (3.36) 

Proof of Theorem \2.1(A We first establish the inequality 

(svp\A(s,I N ,Nx)\) < + E n (A m (In, Nx)) + Eat 




M 


1 


l-M 


da (In, Nx) — 1 
(3.37) 

It has been given in jKS99j . so we only sketch its simple proof. For fixed Sj-i < s < 

Sj with j < M — 1, we have 

A(s, In,Nx ) < \In\~ 1 f da(I]y,Nx ) — G(sj- 1 ) = A(sj,I Nl Nx) + G(sj) — G(sj- 1 ) 
Jo 

= A(sj,I N ,Nx) + jj< A m (In,Nx) + 

Similarly, 


-A(s, I N , Nx) < A m (In, Nx) + —. 


The term 


jj~-y da(lN,Nx) — 1 stems from an analogous estimate for A(sj,Iiy,Nx) 


with s > sm- i- 

Next, we show that 

Eat 


1 


da(Iw, Nx) — 1 


= 0(A~ 


I T I / —v-JV - I -V'W 

l-'ivl Jr / 

Using the notation S(In, Nx) := Jf{i : Nxi € Uv}) we can write 


(3.38) 


Eat 


1 


l-fvl 


da(I]y,Nx) — 1 


= Eat 


S{I n ,Nx) — 1V oin . S(I n ,Nx)~ 1 , 1 

-jTwj- ) - 2Ew -jT^j- + 1 

(3.39) 

Hence, we need to calculate the first and second moment of S(In, Nx). By an easy 
computation, we obtain (using the symmetry of R^) 

^n(S{In,Nx)) = j l/ JV (A r U) N | R^jt i, • • • ,tx)dti... dt N 

Jti,-,t N \ i=1 J 

= jjjj RN(ti/N)dti = |/w| (l + O (J-^j . (3.40) 
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In a similar fashion, using f j2 S 2 (x,y)dxdy = \In\ 2 + 0(\ /at|), we have 


E n (S(I n ,Nx) 2 )= I 
Ju 


N 


2 


ti N (Nti)ti N (Ntj) J Rx(ti,... ,t N )dti ... dt N 

iM =1 


= \InV\i + o 


1 

k n 


+ 0(\I N \), 


which shows together with (13.3911 and (13.4011 


Eat 


\ j n\ 


da(I]\r,Nx) — 1 


= o(4)' 


(3.41) 


Now Jensen’s inequality proves the claim in (13.3811 . We further use the crude bound 


M—l 


En(Am(In,Nx)) < ^ E n (\A(si,lN,Nx)\). 


i —1 


Now, choosing M as the smallest natural number larger than A ]^ 4 , we get with 
(|3.37l) and Corollary 13.41 


Eat (sup\A(s,I N ,Nx)\) < + A^O 


\s£]R 


A 


1/4 

N 


A 


1 / 2 —£ 
■N / 


+ o 


A 


1/2 I ’ 
N / 


which proves Theorem l2.10l for \In\ 1 g(In, Nx). To deduce the result for Nx), 
let for 0 < i < 1 denote 


A := < x 


/ 0 °° da(I N ,Nx) 


I j n| 


-1 


<C A L 


We will assume that N is so large that x € A implies / 0 °° da(lN, Nx) > 0. Then 
da(lN, Nx) — G(s) 


Eat ( sup 


< Eat 1 a(x) sup 


1 


f 0 °° da(I N ,Nx) Jo |-7jv| 


•da(lN, Nx) — G(s) 


+ E n (A c 


It is straightforward to check that x € A implies \In\/ / 0 °° da (In , -Nx) = 1 + 0(A N 1 '), 
where the O term is independent of the specific x. This gives 


Eat ( sup 


f da(lN,Nx) — G(s) 

Jo 


[ — da(I N ,Nx ) - G(s) 

Jo dat | 


+ 0(A n l ) t A (x) sup 


fo da(I N ,Nx) 

l-M 


< Eat ( ^-a(x) sup 

V sGK 

+ E n (A c ) 

= 0(A- 1/4+e ) + 0(A^) + ¥ n (A c ). 

It remains to estimate the probability of A c . The bound (13.4111 gives with Markov’s 
inequality 


PaK^ c ) = o(A“-L). 

Now the theorem follows choosing i = 1/4. 


□ 
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4. Proof of Theorem 12.81 (T1 and (2) 

We need to introduce some of the notation of [KSSV14] , Let us define the 
Mhaskar-Rakhmanov-Saff numbers ay and by via the relations 

[ V , V ' {t) t V , tV ' lt> (4.42) 

Ja v V (by - t)(t ~ ay) Ja v {by - t)(t - ay) 

It is known that for convex, smooth V, ay and by are uniquely determined by 
(14.421) and that these are the endpoints of the support of the equilibrium measure 
Hy. Moreover, it is important for us to see them as functions of V. 

The linear rescaling that maps [—1,1] onto [ay,by] is denoted by 


Ay : 


by - ay by + ay 

v V(s) := --- s +- - -• 


Its inverse is 


A y\t) = 


by + ay 
by — ay by — ay 


-t — 


Hence, [—1,1] C J := A v (J) (cf. (GA)i(l)). Moreover, set 


hy : J x J 


Gy : J 


PV 


hy(t,x):= / (V o \y)"(x + u(t — x)) du 

Jo 

(I/O A y)'(t) - (V o Ay)'(a;) 


t — x 


Gy(x) : = 


1 f 1 hy(t, x) 


£ 


vr i Vl - t 2 


d t, 


pvix) := \l^ Zr ^ 2 Ov(x) 


if |x| < 1, 
else, 


a : (—1,1) —>• R, a(x) := 


1 — x 
1 + x 


1/4 


(4.43) 

(4.44) 

(4.45) 

(4.46) 

(4.47) 

(4.48) 


Note that py is the equilibrium measure of V rescaled such that its support is [—1,1]. 
The actual equilibrium measure py is related to py via 

Pvit) = - — - Py(Ay 1 (t)). 

by — ay 

Proof of Theorem \2.S\ (1) and (2). We will mostly deal with the more involved case 
(2). Let us abbreviate V, f for V — f /N. We will also write || • || instead of || • ||/> 
Combining [ KSSV141 Theorem 1.3] with [ KSSV14 . Proposition 4.1] gives 

h vj-av,l„ „ x 1 (‘('1 , 

—2— Kn.vAXvAOMjO)) = - ^ 


1 

H-cos 

2t r 


N 


g(r,s] 


1 


+ 


1 


a(r) — a(s) 


a(r ) 

+ o 


r — s 


r — s 


(4.49) 


1 ') \a(r) a(s) y 

where g(r, s ) is some function which is not important here. Formula (14.491) holds for 
all r, s € (—1 + 5 ,1 — 6) with arbitrary 5 > 0, where the O term is uniform in r, s 
for fixed 5 and uniform for V — f /N € Xp for some D. In order to use (14.491) for 
the proof of Theorem 12.81 we first have to show that for N large enough and some 
6 > 0 

l-l rir-ii 


A v^Fy L (t/N)) € (-1 + 6,1- 5), for all / with ||/|| < N" (4.50) 

for some p > 0 small. By [ KSSV141, Lemma 2.4], the maps V eA ay, V by, 
defined by (14.421) are Frechet differentiable with (uniformly) bounded derivatives on 
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a neighborhood of V. This lemma was proved in [KSSV14] only for V satisfying 
(GA) 2 but the proof goes through also for V with (GA)i. Hence 

a VJ = a v ( 1 + 0(\\f\\/N)), b VJ = by{ 1 + 0(\\f\\/N)) (4.51) 

and thus 

Xvj(t) = X v mi + 0(\\f\\/N)) (4.52) 

uniformly for t € [—1,1]. For given 0 < r) < 1, assertions (|4.51l) and (14.521) hold 
uniformly for all / with ||/|| < N v , which proves (14.501) . We have thus shown that 
formula (14.491) can be applied. The second summand on the rhs (14.491) is uniformly 
bounded for r, s € (—1 + 6 ,1 — <5) and hence negligible when multiplied by jj. Now, 
we claim that uniformly on (— 1 + 5 ,1 — 5) 

a(r) a(s) l9 , 

^ + ^ = 2 + 0 r-s 2 . 

a(s) a{r) 

Setting z := this claim is equivalent to the relation (z — l) 2 / 'z = 0{\r — s| 2 ). It 
a \ S ) 

is now a straightforward application of Taylor’s formula to show z — 1 = 0(\r — s|). 
Writing 

t 


t:=K 


-l 

V 


N 


s := R 


-l 

V 


N 


we arrive at 


sin 


-KNyjit.s) = 


N7r f^ PV ’ f(r)dr 


+ 0(1/N). 


(4.53) 


(4.54) 


N/j,y(t) ’ ’ N/j,v(t)n(t — s) 

Now, using (14.511) together with the definitions (14.441) . (14.451) . (14.461) and (14.471) . we 
find 

Xvjit)- 1 = Xvitr^l + omm and p VJ (t) = p v (t)(l + <D(\\f\\/N)) (4.55) 

uniformly on [—1,1], Define 

2 f* 

Wj(t) ■= T-—- PvAXvjW) and F v,f(t) := / nv,f(s)ds. 

°v,f a V,f Jav.f 

Furthermore, let gv,f(t ) := Fyj(i) — Fy(t). We conclude 

r X Vj(F) ^ f — g 

/ , Pv,f(r)dr = ( F v + 9vj){t) - (F v + g V j){s) = —— + gyj(t) - gvj(s). 

* 

Using (14.551) and the smoothness of py (see (14.471) 1. it is straightforward to establish 
the relation 

pyj(t) = py(t) + 0(\\f\\/N) 

uniformly on R. It follows that 


gvj(t) -gvj(s) = 0(\\f\\/N)\Fy l {t/N) — Fy i (s/N)\ = 0(\\f\\/N) 


|t — s\ 


N 


where we used in the last step that py is bounded away from 0 on [ay + 5, by — 5]. 
Hence (I4.54|) reduces to 

1 * sin (-kU — s) + C>(||/||/A)|t — s|) 

K N yj(t, s) = - V "V - + 0(1/N) (4.56) 


Np v (t ) 


sin 


Npy(t)ir(t — s) 

(vr(t-s)) 0{\\f\\)\t-s\ 


+ 


Npv(t)Tr(t — s) N 2 py(t)ir(t — s) 


+ 0(1/N), (4.57) 
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where we used the simple inequality |sin(i + s) — sin(f)| < |s|. 

If 1 /|t — s\ = 0(1/N), i.e. > c for some constant c > 0, then Fy 1 (t/N ) — 
FyHs /N) is bounded away from 0 and hence we see that the first term on the 
right-hand side of (14.561) is 0{\/N) which proves the theorem in this case (as the 
sine kernel of t — s is then also 0(1/IV)). 

If \t — s| = o(N), using Taylor’s expansion on Fy 1 * * 5 * * at leads for some v 
between t/N and s/N to 


Fi 


-l 

v 




1 (* ~ s ) _ h F -i (t - 8 ) 2 

li V {Fy\t/N)) N 2 l y M ’ IV 2 ' 


Hence, the second summand in (14.571) is 0(||/||/./V) and we arrive at 


N/j, v (t ) 


k n,vj (t, s) = 


sin(7r(t — s)) 


ir(t - s) - TvA(t, v, N)^X 


+ 0(\\f\\/N), (4.58) 


where 


A(t,u,N) := MFy 1 (tmliFy 1 )"^). 


Note that in (14.58[) the O-term will be of order N 1 for / fixed in case (1) and N v 1 
in case (2). By the simple equality 


1 1 _ -b 

a + b a a{a + b) 


we get 

sin(7r(f — s)) sm(n(t — s)) sin(7r (t — s))nA(t, is, N ) 

ir(t — s) — 7T A(t, IS, N ) Tl(t — s) 7T(f — s)(Tv(t — s) — 7T A(t, v, N ) ) 

sin(7r(t — s)) 1 sm(ir(t — s))A(t,v, N) 

n(t-s) ttN i _ A ( tj l/; jv) (1^) 

= sin(7r(f - s)) f 

7T(/• - ,s) ~ V -Vy ' 


The last equality is due to the boundedness of A(t, v, N) for t € In and our assump¬ 
tion {t — s)/N —> 0. □ 


5. Repulsive Particle Systems - Proof of Theorem 12.81 (3) 

Note that part (3) of Theorem 12.81 for unitary invariant ensembles follows im¬ 
mediately from part (1) and the determinantal relations (|2.17l) . To prove (3) for 
repulsive particle systems, we need to introduce some of the method developed in 

[GV14j to tackle these ensembles. We remark that in comparison to |GV14j . there 
are several new (technical) elements, in particular the truncation to ||/||d < N K and 
the necessity to work with complex-valued processes. Furthermore, aiming at rates 
of convergence requires a separate investigation of the cases of negative-definite h 
and of arbitrary h. 

The first step is to decompose the additional interaction term '^2 i< jh(xi — xj) 
into a linear term and a bivariate term of lower order. This will be done by the 
Hoeffding decomposition w.r.t. a (so far arbitrary) probability measure // on R. 
Setting hfj,(t) := f h(t — s)d^i(s) and for another measure v on R := / Jh(t~ 
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s)dpi(t)dis(s), we may write 

h ( Xi ~ x i) = \Y1 h ( Xi ~ 

i<j i,j 

N 1 

=N^h /J ,(x j ) + 2 \- h ( Xi ~ x J') ~~ Nh v( x i) - Nh^Xj) + /i w ] + C N , (5.59) 
3 = 1 i,j 

where we set Cn '■= — (N/2)h(0) — (N 2 /2)h flfJj . The term N^2^ f =1 h IJi (xj) is of the 
same shape as Q{ x j ), giving rise to the external held := Q + h lt . Our 

aim is to choose p, such that Pjyn and the unitary invariant ensemble Pny^ have 
the same equilibrium measure. To achieve this, the statistic 

Kfi(x) := ^ [h(xi - xj) - Nh ft (xi) - Nh^Xj) + h^] 


should be concentrated under Pn,Vu ■ As U M is a global statistic, we should have 

(5.60) 


lim _ 

AT—>-oo N Z 


^Ny^ld^^x) — 2 (fiw 2/t^ T fa^u), 


where u now is the equilibrium measure to l^. As we will not divide by N 2 and 
thus will be at the scale of fluctuations, the rhs of (|5.60l) should be 0. This leads us 
to the condition is = p, or in other terms, /i should be the equilibrium measure to 
Vp. This recursive problem was solved by a fixed point argument in [GV141 Lemma 
3.1], showing existence of a measure p with the desired property. The uniqueness 
followed later by proving that any fixed point is the limiting measure for Pjf q- From 
now on let p denote the fixed point, set V := and U := The identity 


Pn,q( x ) ~ 


Zn,v 

Zn.vu 


P N , v ( x ) e u W 


with Zn,v,u '■= Znq£ Cn allows to carry many properties from Pn,v over to Pjf q- 
Concentration of U under Pn,v was proved in [GV141 Proposition 4.7] by showing 
that the ratio Zjsry /Z^ypi is bounded in N and bounded away from 0 provided that 
aQ is large enough. More precisely, for any A > 0, there is a constant a(A) < 00 
such that for some 0 < C\ < C 2 < 00 and all oq > a( A) we have 

Ci < E Nye™ < C 2 (5.61) 


for all N. A main ingredient to this is the following concentration of measure result 
for linear statistics (cf. |GV14l Corollary 4.4]), which will be used lateron. 


Proposition 5.1. Let Q be a real analytic external field with Q" > c > 0. Then 
for any Lipschitz function f with third derivative bounded on an open interval D 
containing supp pq, we have for arbitrary e > 0 

E v,Q exp {e ( ^2 f ( x j) ~ N j f(t)dp Q (t)) } < exp { 6 + eCGI/Hoo + II/ (3) ||00)}• 


3 = 1 


Here C is uniform in f, Lip f denotes the Lipschitz constant of f on D and 
is the sup norm on D. 


The key to the local statistics is a linearization method, which transforms the 
bivariate statistic U into random linear statistics. We give an outline for negative- 
definite h, that means h < 0, where h(t) := - 7 = f R e~ zts h(s)ds denotes the Fourier 
transform of h. For such function, —h may be seen as the covariance function 
of a centered stationary Gaussian process (/(i))teR, he. a stochastic process on M 
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whose finite-dimensional distributions are all multivariate Gaussian and such that 
Cov(/(f), f(s )) = —h(t — s ). Then a quick computation verifies that 

1 N 

exp { - 2 Yl h ( Xi ~ = Eex p{5Z ( 5 - 62 ) 

hi i =1 

where the expectation is w.r.t. the probability space underlying the Gaussian pro¬ 
cess. S. Jansen has pointed out to the second author that the linearization (15.621) is 
known in mathematical physics as the Sine-Gordon transformation. Furthermore, 

N 

exp {U(x)} = Eexp{^ f(xj) - N 

3= 1 

holds. The term f( x j)—N f fd/x can now be added to N^jLx V{ x j)i resulting 
in a perturbation of lower order, which does not influence the equilibrium measure. 
The limiting bulk correlations are not altered by the function / either, as can be seen 
from Theorem l2.81 It should be noted that the scaling of the correlation functions is 
independent of /. To summarize, the ensemble Pfc q is an average over determinantal 
ensembles Pj^y—f /N with a small random perturbation of the external field. We will 
show that universality of Py q can be deduced from universality of the invariant 
ensembles P^y-f/N■ Note that the averaging over / results in a weaker rate of 
convergence as uniformity in / has to be shown (cf. Theorem 12.81 ill and (2)). 


fdfj} 


(5.63) 


5.1. Alternative representation of correlation functions and truncation. 

Let us define the generalized invariant ensemble 

1 


P N,Q,f ( x ) 


n 


yM 

N,Q,f l<i<j<N 


. _ x .| 2 p- M Ef=i Q(xj)+J2j = i ffa) 


(5.64) 


where Me N. If M = N, we have Pn,qj = p nq / an d p nq = P nq /> if / = 0. 
Then the k- th correlation function of Pnv at ti, • • • ,t k can be rewritten as 


P N,v(t 1) ■ ■ ■ Pk) 

N\ f 1 


N 


(N - k)\ J^N-k Z N y 


exp 


{~ N ^2 V(xj) + 2 ^2 lo S \ x i - x i\) 


j=k +1 


i<j\ i,j>k 


x exp { - iyy~V(ti, ■ ■ ■ ,t k ) + 2 log \ti — tj|} 

j =1 i<i', i,j<k 

xexp{2 log \ti — Xj\}dxk+i ... dxN 

i<k , j>k 

J\T J 

F (t) exp{2 Z, k)g|L -Xj\} with 


i<k, j>k 


(N — k)\ W Z N y 

k 

F(t) := exp{-IV^V(tj) + 2 ^ log|L -tj\} 

j =i *<i; i,j<k 

Labeling the eigenvalues of the ensemble P$_ k y by x k+ \,..., xy and denoting 

gN 

O := 0%_ k v (t,x) := 2 ^ log|ti-aj|+log [P(t) ], 


i<k , 


(5.65) 
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we arrive at the representation 

R-Nytt i> • • • ;4) = fjfZT /yy®vv-fc,v exp {O}- 


(5.66) 


By analogous steps, we represent the k -th correlation function of Py 


N,Q 


as 


R^ k q{ti ..., tfc) = F j7 ?' u ] Ejv-fc,v ex P (*, s) + O}, (5.67) 

E7v,i/exp{T/(x)) 

where we abbreviated 7/(ti,... ,tk,Xk+i, ■ ■ ■ ,xn) by U(t,x). By (GV14 . Lemma 28] 
we can assume that Xk+i, ■ ■ ■ , xjv G [— L, L] for L large enough. To be precise, the 
lemma shows that for each k we have L, C > 0 such that for all N and for all 
tii ■■■ ? 




^- — - 17 /!' M E 7V-fc,V;L ex P + Ol}| < e 

Ejv,V;Lexp {Z7(z)j 


—CN 


15.681 


where EjJ v . L denotes expectation w.r.t. P^v-h which is the normalized restriction 
of Pffy to [—L,L] n and Ol is the analog of O, obtained by replacing integrations 
over M by integrations over [-L, L\. For later use, let us also state one more inequal¬ 
ity from that lemma, 

N-'R h J Q (t) < exp {CN - ci N[V(t) - c 2 log(l + t 2 )]}, (5.69) 

valid for some constants C, ci, c 2 > 0. It will allow us to restrict the whole ensemble 
(instead of correlation functions) to some compact [— L,L\. 


5.2. Linearization and proof of Theorem 12.81 (3) for negative-definite h. 

Let us give more details on the linearization method for negative-definite h. For 
such h, —h can indeed be seen as the covariance function of a centered stationary 
Gaussian process on M such that (I5.62|) and (15.631) hold. Since the sample paths of 
that process will become a part of the external field, we have to show analyticity. 
This can be done by invoking an explicit representation. Recall that h(t) denotes 
the Fourier transform of h and that we have —h< 0. 

For the representation of /, let (B't )t denote two independent Brownian 

motions and define 

/(f) := (2/V) 1 / 4 / cos(ts)v/— h(s)dBl + (2/tt) 1 / 4 / sin(fs)y — h(s)dBg. 

Jo Jo 

(5.70) 

Here it is convenient to understand the stochastic integral as a Wiener integral, 

/*oo roc roc 

/ 9 {s)dB\ ■=- B l s dg(s) = - B l y/(s)ds. 

Jo Jo Jo 

which exists for g sufficiently smooth and of a certain decay at ±oo (note that by 
the law of the iterated logarithm, \B t \ is almost surely bounded by y/2t log logf). 

Using that Gaussianity of / is equivalent to Gaussianity of all linear combinations 
of the random variables {/(f) : t G R}, it is not hard to check that /(f)teM defined in 
(15.701) forms a Gaussian process on R. Furthermore, it has mean 0 and covariance 
function —h. 

By the assumption on the exponential decay of h, the rhs of (15.701) can be extended 
to a strip {x + iy : x G R, |y| < c} for some c > 0 which implies analyticity of / in 
that strip a.s.. Let 


D := (-L — 5, L + 5) x (-c/2, c/2) 


(5.71) 
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with 6 > 0. Then it also follows from (15.7011 that the extended process {f{w)) w& D is a 
comp lex-valued centered Gaussian process with covariance function (wi) f ( 1 V 2 )) = 

—h(wi — Wf). 

Recall the abbreviation 


S k (t) = det 


sin(7r(tj - tj)) 
ir(ti - tj) 


and set 


(5.72) 


JJ Jl <ij<k 


tj ■■= ( F*)-\tj/N ), C{t) := N~ k niij)- 1 . 

3 = 1 

To prove Theorem 12.81 (3), in view of (15.6811 we have to show 

(iV - fc)! E ^_ fc T; Lexp x ) + °l} ~^N,v-,Lexp{U(x)}S k (t) = O (N~ 1+£ ) 

(5.73) 


for any e > 0 in the prescribed uniformity. Here we used that by (15.6111 , Ejv,v exp { U(x )} 
is bounded and bounded away from 0, which carries over to the truncated setting. 
Note also the slight abuse of notation by not indicating the local scaling of the tj’s 
in Oh- 


Proof of Theorem \2.8\ (3), negative-definite h. Let h be negative-definite. Further¬ 
more, let / denote a centered Gaussian process with covariance function — h and 
define /:=/—/ fdfi. Then we have 


TV 1 

C(t) ' wN 


=C(t) 


(N — k)\ 
N\ 

(N - k)\ 


E N-k,v-,L ex P { u {i, x) + O l ] 


N 


E 


E N—k,V\L ex P {Y1 /((*> X )i) + °4 


3 = 1 


N 


=E 


e n,v-,l exp (E/P;)}C(t)<K/; L (4 1 ■ ■ ■ j 4 ) 

3 =1 


where we used the identity 

P k ,* n N\/(N — k)\ 

K N,VJ-dtl, ,4) - --E 

e n,vj-l^= 


N 


E N—k,V;L ex P {Y /((*> X )i) + ° L } 
3 = 1 


which can be obtained analogously to (15.6711 . Thus the first summand of (15.731) 
equals 


E 


N 

E N y. L exp{Yf( x j)} {c{t)R k N yj. L {ii,... ,t k ) - Sk 


3 =1 


(5.74) 


To apply Theorem l2.8l f2h we will replace the integration over all / by an integration 
over / with ||/||r> < N 71 , where D is the complex domain defined in (15.711) . More 
precisely, we will show that 


N 


El 


d>Nv} 


e n,v-,l exp 


3= 1 


X (C(t)R%yj.y(t\, . . . ,4) — Sk 


= 0(e~ cN2v ) 


(5.75) 
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for some c > 0. This will be done by applying Holder’s inequality to separate the 
expectations of 1 {||/|| d >a^}> of Eyyy exp { Xy=i f( x j)} an d of 

C(t)R%yj. L (t 1, ■ ■ ■ , tk) — S k (t). 

Let us reconsider the complex-valued process (f{t))teD, which is the extension of 
the Gaussian process / on [—L,L\ with covariance function —h. It follows from 
(I5.70P that real and imaginary parts of / on D are (real-valued) centered Gaussian 
processes. Their covariance functions are readily computed as 

E(Re f(wi) Re /(u^)) = --(Re h(w\ — W 2 ) + Re h{w\ — W 2 )), 

E(Im f(wi) Im /(u^)) = -(Re h(w 1 — W 2 ) — Re h{w\ — W 2 )), w\,W 2 € D, 

giving the variances 

E(Re/(w)) 2 = --(Re h( 0) + Re h(2ilmw)), 

E(Im/(u>)) 2 = ^(Re h( 0) — Re h(2ilmw)), w € D. 

Now Borell’s inequality (see e.g. ]AT07l Theorem 2.1.1]) states that the supremum 
|| X || k of a continuous centered Gaussian process X t over a compact K has sub- 
Gaussian tails, more precisely it is dominated by a Gaussian random variable with a 
certain expectation and variance a := sup teK EX] 2 . Since the sum of sub-Gaussian 
variables is sub-Gaussian as well, we see with sup, i( , sdI/ + MI < sup w&D |Re/+(u;)| + 
su ProeB|I m f + ( w )\ that sup wgD |/ + (t(;)| is also sub-Gaussian. The same reasoning 
leads to the conclusion that sup wgD |/(u;)| is sub-Gaussian, giving 

P{\\f\\ D >N*} = 0(e- cN * n ) ( 5 . 76 ) 


for some c > 0. 

Next, we provide a bound for Ejv,v ; l exp { f( x j)}- Proposition 15.11 also 

holds for the ensemble truncated to [— L, L\ with an exponentially small error which 
we will neglect. However, it is crucial that in the truncated case the Lipschitz 
constant is taken over [-L, L] instead of the whole real line. Thus we get 

E n,v-,l exp { ^2 /( x i)} ^ exp { + C (WfWl ~ L M + ll/ (3) ll[-i-,i])} ( 5 - 77 ) 

3 = 1 


for some C. f is again a centered stationary Gaussian process with covariance 
function —h" and thus, by similar arguments as above, sub-Gaussianity of Lip/ and 
analogously also sub-Gaussianity of ||/^ II[-l,l] follow. Hence for some A > 1 (close 
to 1, coming from Holder’s inequality) there is a constant C such that for all N 


E 


E n,v-,l exp 

3 = 1 


<C, 


(5.78) 


if olq (and hence ay) is large enough. 

It is important to note that, as the processes are stationary on M, the variances 
of Lip/, ||/1 | { _ LM and ||/ (3) ||[_ LiL] and therefore the required are independent of 
the truncation threshold L and k. 

Now we will estimate C(t)R k N v f- L (h, ■ ■ ■ ,t k ) — S k (t). If f* = tj for some i / j. 
then this quantity is 0, hence we will only consider t with distinct elements. For 
such t, (12.171) and Pn,vJ;l(. x ) > 0 for any x € [—L,L\ N with distinct components, 
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imply that {KN,vj-,L{Uitj))\<i,j<k =: A is a positive definite matrix and can hence 
be written as A = B 2 for some matrix B. Using Hadamard’s inequality then gives 

k k k 

det A = (det B ) 2 < JJ \B 

i=li=l j=l 


which is 

k k 

■ > 4 ) < K N ,VJ-L(ij,tj) = RN,V,f;L(ij)- (5.79) 
i=i i=i 

Let us employ the representation as inverse Christoffel function (see e.g. [TotOOp 
d i ^ e -NV+f 

RN ’ v ’f {t) = A N (e- NV +f,tV 

X N (e- NV+f ,t) := inf / 

Pzv —1 (*)=1 J 


where the infimum is taken over all polynomials Pn-i of at most degree N — 1 with 
P/v_i(t) = 1. This representation immediately implies the bound 

< «W(4> 2|I/I11 -^- 

To bound R l NV . L , we can now use the uniform convergence stated in Theorem 12.81 
(2) together with the boundedness of the sine kernel, giving 

N~ l R} N V j. L {tj) < Ce 2 ^^~ L P and hence 

C(t)R k N yj. L {iu ... , 4 ) < c’e 2k m-w, (5.80) 

where C,C' do not depend on /. As ||/||[_l,l] is sub-Gaussian, we have for some 
A' > 1 and some C" that for all N 

E|C7(t)i2^ )V|/;i (ti,... , 4 ) - S k (t)\ y < C". (5.81) 

Combining (15.761) . (15.781) and (15.811) . it follows that (15.751) is of order 0(e~ cN2ri ). 
Thus Theorem 12.81 (2) yields 


N 


E 


Ejv,F;L exp 


V,f;L 


ik)~ Skit)) =0(N~ 1+£ ), 


3 =1 


for any e > 0, given that olq is large enough. 


□ 


5.3. Proof of Theorem 12.81 (3) for general h. 

A general h may be decomposed into positive-definite functions as follows. Let 
± denote positive and non-positive part and write h = ( h) + — (L)_. Then h = 

h + — h~, where h± := [h)± and furthermore, are positive-definite. To use 
for our linearization method, we need these functions to be real-analytic. The 
real-analyticity is somewhat surprisingly equivalent to the exponential decay of h 
at infinity. On the one hand, exponential decay of h allows to extend h± to the 
complex plane via Fourier inversion, thereby showing real-analyticity. On the other 
hand, [LS521 Theorem 2] tells us that any real-analytic positive-definite function has 
a Fourier transform of exponential decay. 
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Instead of —h, we will use —h z := zh + + h , z > 0, as a covariance function and 
use complex analysis to show the desired. To this end define for complex z € C 

N 

Uz{x) :=|( Y2 h+ ( x i - x i) - [fyU**) + h $( x i) - Kit] ) ( 5 - 82 ) 

i,j= 1 

| N 

+n ( h ( Xi ~ x i) “ [fyi ( x *) + h n ( x i) - v] ) • ( 5 - 83 ) 

i,j =1 

Again, we have to show 

(^ C ^)! E jv-fc,v ; L ex P x) + 0 L } -Ej Vi y ;L exp{W z (a;)} < S , fc (f) = o(l) (5.84) 

for z = —1. The proof for negative-definite h implies (15.841) for z > 0. This is 
basically enough, as Vitali’s theorem implies, which we state for the convenience of 
the reader (cf. [Tit39l 5.21]): Let g n (z) be a sequence of analytic functions on a 
domain U C C with \g n (z)\ < M for all n and all z € U. Assume that lim rwoo g n (z) 
exists for a set of z having a limit point in U. Then lim^^oo g n (z ) exists for all z in 
the interior of U and the limit is an analytic function in z. 

The set containing a limit point will be chosen as a small interval (0, 5) for some 
<5 > 0. We remark in passing that 5 can be arbitrarily small and as a consequence 
h + has no influence on the necessary size of the convexity constant cxq. 

To transfer the required uniformity in the tj from the case z>0toz = — lisa 
technical issue as taking absolute values and suprema would destroy the analyticity 
in z, which is necessary for the application of Vitali’s theorem. Therefore, we will use 
the following characterization of uniform convergence in terms of sequences: A se¬ 
quence of continuous real-valued functions (g n )m defined on M z , converges uniformly 
on the sequence of compact sets (B n ) n , B n C M. 1 towards a continuous function g if 
and only if for all sequences (n m ) m C N with limm^oo n m = oo and all sequences 
{im)m with t m € B nm we have lin gn rn (trn) g{^m) — 0. 

We will take B n := 1^- Let (N m ) m C N be a sequence going to infinity and 
(' tt m ))m be a sequence with t< m ) € INm- Let us define W m : C —>• C as 

W m (z) :=C (f( m ))IEjy™_^y;£ exp , x) + Ol} 

- E Nm y-,L exp {U z (x)}S k (t im )). (5.85) 

Note that we suppressed some of the m-dependencies. Clearly, (W rn ) m is a se¬ 
quence of analytic functions. 

Remark 5.2. Vitali’s theorem allows to deduce convergence in a region of the complex 
plane from convergence in another region. As rates of convergence might well be 
different according to the region one is looking at, it is clear that we cannot transfer 
the rate 0(IV~ 1+£ ), valid for z > 0, to z = —1 with the same technique. This seems 
to be rather a technical issue, we in fact believe that the correct rate should also be 
(at least) 0{N~ l+e ). 

Proof of Theorem 1 2. fll (3), general h. We will often drop the dependence on m in 
the following. For positive z, we can apply the linearization procedure as described 
above, as then h z is a positive-definite function. Thus we have W m {z ) = o(l) for any 
z € (0, 5) for some 5 provided ccq is large enough, where 5 is chosen so small that 
the lower bound on ctQ does not depend on <5. Hence, to apply Vitali’s theorem, we 
need to show boundedness uniform in m and 2 from a domain containing —1 and 
(0,5). This domain may in fact be chosen as the halfplane {z € C : Rez < 5}, 
which can be seen as follows. Bounding W m termwise, (15.820 shows that Im z only 









EMPIRICAL SPACINGS OF UNFOLDED EIGENVALUES 


29 


gives a phase which vanishes by taking absolute values, hence we can concentrate 
on real z. For z < 0, (15.821) is non-positive, since it is minus 1/2 times the variance 
of a Gaussian random variable (cf. (15.631) 1. implying that in this case the influence 
of (15.821) in bounding W m will vanish as well, due to taking absolute values. It is 
thus sufficient to consider z > 0. In this case we get as above that W m (z) equals 

N 


E 


Ejvy ; i exp { fz(xj)} (c(t)R%yj z . L (t 1 ,..., 4) - s k 


(5.86) 


3 = 1 


where f z is a centered Gaussian process with covariance function — h z . Now, the 
bounds (15.731) and (15.801) can be used again to show uniform boundedness in m. To 
check that these bounds are uniform in z € (0, 5) for 5 > 0 small, is straightforward 
and left to the reader. The theorem is proved. □ 


6. Proofs of remaining statements 

Proof of Theorem 1 2. ‘A Let an interval In C [0,1V] be given such that 

lim inf n^oo dist(/jv, {0, N})/N > 0. We will apply Theorem 12. 101 with Pn being the 

distribution of the unfolded ensembles Pn,VJ an d Pn q, respectively. We start with 

the former case. Let Fy ^ be a function with the following properties: 

(1) Fy 4 R A J strictly monotonically increasing, continuously differentiable, 

(2) 4" 1] (M) = J, 

(3) Fy ^(t) = Py 1 ^) for t € U with U D liminfjv->oo In/N open. 

Now, we define 

N 

P N {x) := HFl- 1] \x j )PNy,f(Fl- 1] (x 1 ),...,Fl- 1 \x N )). 

3 =1 

By properties (1) and (2), Pn is a probability measure on M. N . By property (3), we 

have F [ ~ 1] (t) = {nv{P\ 4 ( 4)) -1 f° r all t suc h that Nt € In for N large enough. 
Then the correlation function /?(/ is (for N large enough) on /(/ given by 

k 

R k N (h/N,... ,t k /N) = n MFv\h/N))- l R k Nyj {F- l {ti/N ),..., F-\t k /N)). 

3 = 1 

It remains to check conditions (1) and (2) in Theorem 12. 101 For condition (1), we 
may use inequality (15.801) . giving 

k k 

Y[n v (i j )~ 1 R k N y f (ii,...,i k ) < 
j = 1 i =1 

which is valid also for ensembles with non-compact support. Here we use again the 
abbreviation tj := Fy 1 (tj/N). Positivity of fiy on (ay,by) yields condition (1) for 
the unitary invariant ensembles. 

For condition (2), another application of Hadamard’s inequality may be used, 
which we cite from [AGZ101 Lemma 4.3.2], For kernels K\,K 2 , defined on some 
locally compact A x A and all t € A k , the following inequality holds: 

det/Ad (L, det(A 2 (L, 4'))l<i, j<k\ < fc fc/2+1 1|ATi - AT 2 || • max(||A' 1 1|, ||A' 2 ||) fc 

(6.87) 

Here ||A'|| := sup t sey ^\K(t, s)| denotes the sup-norm of a kernel K on Ax A. Choos¬ 
ing Ad(M) := (lV//u(A^ 1 (t/lV)))- 1 iL Nyi/ (A- 1 (t/Ar),A- 1 ( S /N)) and K 2 as the 
sine kernel, we find with A := In that by Theorem 12.101 (1) \\K\ — K 2 \\ = 0(1/N). 
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Moreover, the boundedness of the sine kernel and the uniform convergence of K\ 
imply that max(||iFi||, ||i^ 2 ||) < C for some C > 1. This proves condition (2) of 
Theorem 12.101 for Pn,vj- 

For the repulsive particle systems P^q, we will first truncate the ensemble to a 
compact [—L,L\ n and apply Theorem 12.101 to the truncated ensemble. This will 
be convenient as the truncation threshold L in (I5.68p depends on the order of the 
correlation function k and in the subsequent linearization several constants depend 
on L, which would complicate checking conditions (1) and (2) of Theorem 12.101 
significantly. 

For any L > 0 we have the crude bound 




< N 


s: 


[ n-L,L]")c(x) sup 

sSK 


J-J°dv(I N ,x)-G(s) 


<NP^ q (([-L,L] n Y) 


Now, 115.691) tells us that choosing L large enough, the last expression is bounded by 
exp(—IV), hence negligible. This also shows that replacing the normalizing constant 
q by its truncated variant q. l , obtained by replacing the integrations in 
Z^j q by integrations over [—L,L\, only results in an error of at most exp(— N). 
Summarizing, we get 


e n,q I sup 
\seR 


l-Twl 


= e n,q-,l I sup 

vSEM 


f da(I N ,x) — G(s) 

Jo 

i r s 

da(I]\r, x) — G(s ) 


\ j n\ 


+ 0(e 


-N\ 


where E %q. l denotes expectation w.r.t. the truncated ensemble P^q. l - Note that 
both ensembles have the same limiting measure /jlq and therefore, there is no change 
in the unfolding of the particles. The case of a instead of |/at|“ 1 ct is completely 
analogous. 

The same procedure as above can now be applied to fit Py q. l into the frame¬ 
work of Theorem 12.101 Condition (1) then follows using the linearization procedure 
introduced in Section 5. Indeed, an inequality of the form 


■ ■ ■ i^k) — P 


for some C, uniformly in t\,... ,tk, is equivalent to 


Ejv_fc,y ; L ex P x ) + Ol } < C' 


( 6 . 88 ) 


uniformly in t±,...,tk for small positive z. By linearization, the l.h.s. of ([6.881) is 
equal to 


E 


N 

Eat,v;l exp { ^ 4 ) 

3 = 1 


to which now (15.80|) may be applied for almost all /. The sub-Gaussianity of ||/||oo 
gives the desired bound. 

For condition (2), we can use the same arguments and (16.871) . □ 
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Proof of Corollary \2-4\ Let In be a sequence of proper compact sub-intervals of 
[0,1V] with In/N —» [0,1], We will use the simple inequality 


E N sup 


1 


N- 1 


da([ 0, N], x) — G(s) 


'o 


< E n sup 


1 


N-1J o 
The last term is 

1 


da (In, x) — G(s) 


+ E n 


1 


N- 1 


AT- 1 


S j : Xj G 


_i /[0,JV] \ Jjv 


IV 


du([0, IV] \ In,x), (6.89) 


- 1 


where F = Fy or F = Fq, respectively. F 1 (([0, A^] \ In)/N) consists of the two 
intervals (—oo ,ay + d N ) and (by — en,oo) for two sequences of positive numbers 
£N,d N , both converging to 0. Let us exemplarily deal with (by — en, oo). We have 


1 1 f°° 

— —-E N#{j : Xj € (by - £at, oo)} = ——- / R$f(t)dt 
iV - 1 iV - I Jbv-CN 


(6.90) 


For v j, [KSSV141 Theorem 1.5] provides the following asymptotics. For a certain 
c > 0 and all t G J with t > 1 + c _1 fV“ 2 / 3 , we have 

-ft] 


2ir(byj — ayj)N 


0 - N VV,f{t) 


t 2 - 1 


+ o 


N(t - l) 5 /2 


+ o 


where 


rt 

r]y{t) ■= J Vs 2 — lGy(s)ds 


, (6.91) 


(6.92) 


and Gy has been defined in (I4.46|) . Similarly to (|4.51|) . (14.521) and (14.551) . we can 
show rjyj = rjy( 1 + G(||/||/1V)). As V is strictly increasing and Hindoo V(t) = oo, 
we have for all t large enough and s G [—1,1] that (V o A y)'(t) — (Vo Ay)'(s) > d 


for some d > 0. Thus (cf. (14.45 1) 1 hy(t,s ) > 


t — s 


and (cf. (14.46 j) ) 


Gy(t) > 


n J -1 


7T 7-1 (t - S )Vl - 


zds = 




where the last equality is due to t > 1. Hence r]y(t) > dt and we conclude that for 
some d' > 0 


f 

J bv 


~ °(1)- (6.93) 

lb v +c"N ~ 2 /3 

With the linearization technique, this bound can be transfered to the ensemble P]f q. 
We note in passing that due to the mixing over /, for Pjf q the bound for the first 
correlation function is somewhat worse than (16.911) . see [KV15j for details. 

For the edge regime, i.e. by — en < t < by + c"fV~ 2 / 3 , the following asymptotics 
are given in [KSSV141 Theorem 1.5], valid for 1 — c < t < 1 + clV -2 / 5 (with the same 
c as above), 


Rn. 


1/3 

,v,/(V,/(/)) = _ a v ' f f KAi (jv,fN 2/3 (t - 1), lyjN 2 ^ 3 (t - 1)) (1 + r(t)), 











































EMPIRICAL SPACINGS OF UNFOLDED EIGENVALUES 


32 


where 7 yj '■= 2 1 // 3 GV,/(1) 2 / 3 , 

(0(l--t) + 0(N- 2 / 3 ), iff<l 

7 ^ “ je>(lV(f - I) 5 / 2 ) + 0{N~ 2 / 3 ), if t > 1 

is uniform in / as above and as before, we can neglect the /-dependence. Here K ai 
is the Airy kernel, 

jy Ai(t) Ai'(s) — Ai'(t) Ai(s) 

t — s 

On the diagonal, the Airy kernel is given as K\[(t, i) = Ai/f ) 2 — tAi(f) 2 . We are 
therefore left to estimate 

^' k m (j v N 2 / 3 (\y\t) - l),wN 2 / 3 (\y\t) - 1)) . 

Let us first consider the case t < by, i.e. Ay 1 (t) < 1. Setting ((t) := |f 3 / 2 , [AS641 
10.4.60, 10.4.62] provide the following asymptotics, valid as t > 0, 

Ai(-t) = vr -VH-V* (sin (C(t) + J + 0(C(i) _1 )) , (6-94) 

Ai'(-t) = - 7 T _ 1 / 2 f 1/4 (cos (C(t) + ^) + ^(C^) -1 )) > (6.95) 

from which we conclude 

Ai(-t ) 2 = Oit- 1 / 2 ), Ai'(-t ) 2 = 0{t 1 / 2 ). 

Thus we get uniformly for by — sn < t < 1 

iV-V3^ A . ( 7v/iV 2/3 (A -l W _ 1), 7 vA t2/3 (A -1 W _ !)) = 0{\y\t) - 1) = 0(1). 

For 1 < t < 1 + cN~ 2 / 5 , the formulae [AS641 10.4.59, 10.4.61] provide 
Ai(t ) 2 = 0(t-V 2 e-l t3/2 ), Ai'(f ) 2 = 


which yields uniformly for by <t < by + c"N 2// ' 5 

R] V,v(A) = °(1) 

in an analogous way. This gives 


/ 


6y+c"AT- 2 / 5 


b V ~£N 


RN,V,f{t)dt 


0 ( 1 ). 


(6.96) 


This estimate can by the now familiar procedure be extended to the ensemble q 
Altogether we have combining (16.901) . (16.931) and (16.961) 


jj—jE N #{j ■ Xj € {by - £iv,oo)} = o(l). 


Returning to ([6.891) , we wish to apply Theorem 12.101 for an interval In which ex¬ 
hausts [0, N] for N —» 00 . It suffices to show uniform convergence of the correlation 
functions towards the sine kernel, as conditions (1) and (2) of Theorem 12.101 will 
then follow exactly as in (15.801) and (16.871) . Hence our task is to extend the proof of 
Theorem 12.81 to regions close to the edges. 

Formula (14.491) was valid for r,s € (—1 + 4,1 — 6) for some arbitrary but fixed 
5 > 0. As now I A covers some of the left and right edge regions, we also have 
to consider correlations between particles from different regions. Here the general 
[KSSV141 Theorem 1.3] is useful which gives 


K N,V,f{t, S ) 


ki(t)k 2 (s) - k 2 (t)ki(s ) 


+ 0(N~ 1 ), 


t — s 






























EMPIRICAL SPACINGS OF UNFOLDED EIGENVALUES 


33 


where k\, are bounded functions and the O term is uniform for t, s from bounded 
subsets of J. This means that Kjyyjit, s) is bounded if t — s is bounded away 
from 0, which includes the case of one particle at the left and the other at the 
right edge. Dividing by N, we see that Theorem 12.81 (1) and (2) hold in this case. 
The correlations between bulk particles and edge particles are more subtle, as these 
regions are adjacent and a transition from sine kernel to Airy kernel statistics occurs. 
Without loss of generality we will consider this transition at the right edge only. To 
state the analogue of (14.491) at the edge, we need some more notation. Let us remark 
that we restrict ourselves to the case t < by, the void region not being of interest 
here thanks to (|6.93l) and (16.961) . 

Recalling (14.471) . define 

: M —>• M, £y(t) := 27r py(s)ds 

and for some 0 < <5 < 1 

fo \ 2/3 

}n,v '■ [1 — <5,1] -A M, fN,v(t)'■=-N 2 / 3 (-£y(t)J , (6.97) 

/o \ 1/6 

dy : [1 - 6, 1] -A- M, dy(t) := • (6-98) 

Abbreviating fjy := fyy, [ K3SV14 , Proposition 4.1] states that for some 5 > 0 
uniformly for r, s € [1 — 5,1] 


byj - a V j 


K N,v,f{^V,fir),\yj{s )) = K Ai (f N (r),f N (s )) 


fN(r) ~ fN(s) 


+ 


( Ai(f N (r)) Ai'(f N (s)) 




dy(s) 


+ 


r — s 

Ai (/N(s)) Ai'(/jv(r))\ dy(r) — dy(s) 


d v (r) 


r — s 


+ 0(1/N). 


Note that this formula covers a part of the bulk region. By (|6.94l) and (|6.95l) . 


KAi(fN(r),f N {s))(f N {r) - f N (s )) 

= Ai(f N (r)) Ai'(/jv(s)) - Ai ’(f N (r)) Ai(f N (s)) 


(6.99) 


1 

7T L 

— cos 

1 


cos 


N 


W)- 


7T 


COS 


N 


>)) + 


7T 




(I - ))-t) c °s(C(/iv(s)) +j)j +0 (C(|/iv(f)|) 1 + C(I/jv(s)|) X ) 


= -sm(((f N (r))-((f N (s))) + 0(({\f N (r)\) 1 + C(|/iv(s)|) X ) 


7T 

1 


= — sin ( Nit 
7r 


u) +0 


rr 

' / py{u)dt 

J S 

Hence we have for r,s<l, r/s, 
byj ~ a v,f 


gy(r) ^^(g) 
N 


-i 


- K N,vj(^vj(r),\yj(s)) = 


1 sin (ALr fj pyj(u)du ) 


7r 


r — s 


+ 


( Ai(f N (r))Ai'(f N (s)) + M{f N (s)) Ai'(f N (r))\ dy(r) - dy(s) 


dy(s) 


dy(r ) 


r — s 


( 6 . 100 ) 


+ o( + 0(l/N). 

\ N (r - s ) J 


We have thus recovered the leading term of (14.491) . However, the 0-term involving 
^(r)^ 1 will only be small for r and s being not too close to each other. We will 
therefore first consider the case of |r — s| > N~ 2+p for some small p > 0. From 
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(I6.94|) , (16.971) , (16.981) and the boundedness of the derivative of dy it follows that for 
1 — 6 < r, s < 1 — e' n with e' N > 0 converging to 0 slowly enough, we have 


h v,f - a VJ 

2N 


K N,v,f{^v,f{r), A v,f(s)) 


sin (Nn j; p VJ (u)du) 

vrlV(r - s) 1 ’ 


for some t > 0, uniformly for r, s with |r — s| > N~~ 2+p . Now we can proceed as in the 
proof of Theorem 12.81 (1) and (2). We get with properly chosen e A > 0 converging 
to 0, that uniformly in t, s € [(1 — S)N, (1 — £n)N] 


1 


k n,v.j 




sin(7r(t — s)) 

7T (t — s ) 


+ O ( N ~ L ) . 

( 6 . 101 ) 


Here we already replaced Fyj and pyj by their counterparts Fy and py, which has 
been justified in the proof of Theorem 12.81 

For |r — s| < N~ 2+p , we may use Taylor’s expansion in s at the point r in (16.991) 
together with Ai "(t) = t Ai(t) to obtain 

KAi{fN(r), f N (s))(f N (r) - f N (s )) 

= /v( r )/A r (p) Ai(/jv(r)) 2 (s - r) - f' N (r) Ai \f N {r)) 2 (s - r) (6.102) 

+ \ Ai (/v(p)) [f' N ( v ) 2 A H/v(u)) + f N (v)f'f,(v) Ai (f N (v)) 

+ /v(u)/^(u) 2 Ai\f N (is))] (s - r)(s - u) (6.103) 

+ i Ai'(/^(r)) [/jv(A)/k(A) Ai(/jv(A)) + /&(A) Ai'(/jv(A))] (a - r)(a - A) (6.104) 

for certain u, A between s and r. Using ) = Ai^f) 2 — tAi(f) 2 , we see that 

(J6.102D equals 

f' N (r)K Ai (f N (r ), f N (r))(r - s ). 

^From (16.941) . (16.951) . (I6.97|) and \r — s| < 1W 2+P , we find 

bV ’ f 2N VJ KN ’ Vj{Xvj{r )’ A ^ (s)) = ^ A Ai(/v(r),/v(r)) + 0(1V 1+P ) 

with the 0 term being uniform in r, s. Using (16.941) and (16.951) again, we can for 
t -A — oo derive 

= i(-i)V 2 (i + 0((-t)- 3 / 2 )). 

7r 

With these asymptotics and (16.971) we see that 

^P-K Ai (f N (r), f N (r )) = py (r) + 0(Jr‘) 

uniformly for 1 — A<r<l — e)y with e)y —>• oo slowly enough, which establishes 
(16.1011) also close to the diagonal. As written above, this suffices to apply Theorem 
12.101 to finish the proof. The transfer to correlation functions is made with (12.171) 
and to P A q with the linearization method. This proves Corollary 12.41 □ 


Proof of Corollary \2.1R The corollary follows from Theorem 12.101 by setting 


Pn(x ) := py{a) N Pn,vj ( « + 


Xl 


py{a)’ 


! a 


XN 

p v (a) 
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The case of repulsive particles is analogous. By the two-fold application of Hadamard’s 
inequality (cf. (15.7911 and (16.871) 1 and the linearization method, it suffices to show 

I f t s \ sin(7r (t — s)) fl + |t| + |s| \ 

Nfj,v(a) N,V,f N^y(a)' aJr Nii v {a)) ir(t — s) \ N )' 

This is precisely the statement of [KSSV14 , Theorem 1.8], which can also easily be 
deduced from (14.491) . Here we used again that the /-dependence in the scaling can 
be neglected. □ 


Proof of Corollarv \ 2. (A A careful reading of the proof of Lemma 13.21 shows that we 
have for any s, L and any e > 0 



1 

l^ivl 


dK^a(lN,x) — G(s) 


L 




S “ C » -O 


(*-!)! 



+ 


gk—lfjkj^k/2+1 

(k-iy.N 1 -^ 


+ E(s, L) + E(s, L + 1), 


where E(s, k) has been defined in (13.291) . Choosing L E N such that ydogf/jvl 
and L = o(log|/jv|), we get similarly as in the proof of Lemma 13.21 
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for any e' > 0, where the O term is uniform in 0 < s = 0(ylogj/]v|). Thus it remains 
to see that taking the supremum over [0,0(yTog|ijv|)) is sufficient. It follows from 
the sub-Gaussian tails of G that 1 — G(0(-y / log|/jv|)) = 0(\In\~ 1 ) (cf. (13.361) 1. From 
(13.401) we know that the expected total mass of |/ at |~ 1 o '(/ at , x) is 1 + 0{\In\~ 1+£ ) 
for any e > 0 and hence the uniform approximation on [0, 0(yTog|//v|)) gives that 
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GHVloglbvl) \ T n\ 


dE N a(lN,x) = 0(|ijv| 


-1+eL 


for any e' > 0. This proves the corollary. 
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